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ON THE SOLUTIONS OF ORDINARY LINEAR HOMOGENEOUS 
DIFFERENTIAL EQUATIONS OF THE THIRD ORDER 

Bv George D. Bikkhoff 

Introduction. Given a linear differential equation of the third order 

(1) y"> + py" + qy' + ry = 0, (alasSJ), 

in which the coefficients^?, q, r are real functions of x, continuous together 
with their derivatives of all orders, but otherwise unrestricted ; what is the 
general character of the solutions and how do they depend on these coeffi- 
cients? This rather interesting question lies almost wholly untreated, and it is 
the aim of the present paper to consider it. 

The analogous question in the case of an equation of the first order 

is at once answered. By means of the explicit form of the general solution 

y = c e~ Spdx , (c, a constant), 

we see first that y cannot vanish, unless it does so identically. Also any func- 
tion y x which does not vanish is the solution of such an equation, namely 

y'- y ly= 0. 

Thus the general characterization of the solutions has been obtained. The 
exact nature of the dependence of the solutions upon the coefficient j) is exhib- 
ited by the explicit formula. 

If the equation be of the second order, 

y" + py' + qy = o, 

the genera] solution has the form 

c iyi + c 2 ?/ 2 , (c lt c 2 , constants), 

where y x and y 2 are any pair of linearly independent particular solutions. It 

(103) 
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is an immediate consequence of tho fundamental existence theorem for ordinary 
differential equations that 

W=y 1 yi-y i y' l 

cannot vanish. For suppose that this expression did vanish for x = £. In 
this event we could determine a and /8, not both zero, so that the two equa- 
tions 

«yi(f) + /8y,(f) = 0, ayi(f) + /8yi(f) = 

would hold. The solution 

y = ay x + /3y 2 

and its first derivative would then vanish at x = ij . But, by the existence 
theorem referred to, there is but one such solution, namely y = 0. Accord- 
ingly y x and y. z would be linearly dependent, contrary to our hypothesis. 
Thus W cannot vanish. Furthermore if T^andy^are any functions for which 
W ^ 0, these functions will form a pair of linearly independent solutions of 
the linear differential equation of the second order 



= 0, 



the coefficient of y" being precisely —W. This coefficient can be made unity 
by division of both members of the equation by — W. 

The analytic condition W jt 0, which has been shown to characterize y\ 
and y. 3 , tells us that the ratio of y x to y. 2 increases or decreases throughout the 
interval considered. Let now y x and y. t be regarded as the coordinates 
of a point in the projective line. As x increases through the interval, the 
point will describe the line in one sense continually. The elementary Sepa- 
ration Theorem concerning the roots of y, = and y t — is an evident con- 
sequence of this fact : the roots of y x = and y 2 = will separate each other. 
The above results give a characterization of the solutions, and thus answer 
the first part of our question for a linear differential equation of the second 
order. 

As the most important of the theorems which deal with the dependence 
of solutions upon the coefficients we cite the following Comparison Theorem : 
suppose tAvo equations 

y" + liV = ° and y" + WJ = 0, q x > q 2 , 
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to be given ; then there will be at least one zero of any solution of the first 
equation between two zeros of any solution of the second, or less precisely, 
an increase of q in the equation 

if + qy = 

brings the zeros of the solutions closer together.* 

Very analogous theorems exist in the case of the equation of the third 
order, and it is with these that we are concerned here. 
1. Geometrical Interpretation. If one writes 



W 



a necessary and sufficient condition that i/ lt y 2 > I/s form three linearly inde- 
pendent solutions of some linear differential equation of the third order (1) is 
seen to be W V 0, just as the analogous fact was seen in the case of the second 
order equation. An explicit formula for W is obtained by noting that one has 
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by (1), so that Wis given by the formula 

(2) W=ce-s» dx , ( C * 0). 

*Take for example the equations y" -f-X,i/= and y" + \^y = 0, where Xi and X 2 are 
positive constants, Xj greater than X2. The interval between zeros is fixed in both cases, 
tt/V^Xi in the first and ir/^M in the second case. 

The theorem is more general than at first might appear since any equation of the second 
order may be given the normal form, 

y" + qy = 0, 
by a simple transformation. 

The reader is referred for the proof of this and similar results to articles by Professor 
BScher in the Bulletin of the American Mathematical Society, vol. i, pp. 295-313, pp. 3G5-376, 
1897-1898, as well as to the original paper of Sturm, Journal de Hathematique, vol. 1, pp. 100- 
180, 1830. 
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To a set of linearly independent solutions y ly y 2 , Its of a given equation 
(1), whose general solution is then 

<Wi + <W2 + c 3 y s , 
there will correspond an integral curve C defined by the equations 

Wi = yi (x) , w 2 = y 2 (x) , w 3 = y 3 (x) , 

in which w t , w 2 , to 3 are the homogeneous coordinates of a point in the pro- 
jective plane referred to a certain triangle of reference. It is clear that any 
other integral curve D obtained by making a different choice of the three lin- 
early independent solutions will be projectively equivalent to C, since the new 
solutions can be expressed linearly in terms of y ly y 2 , ?/»• 
The equation of a tangent line to O is 
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= 0. 



The homogeneous coordinates of this line are therefore 



(3) 



Z\ = 



_ y^A - yiy'i y _ ?hyi - y^y* _ _ yiy'* - y%>A 
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the divisor W being so chosen that z x , z. if z 3 are solutions of the adjoint 
equation to (1) : 



(4) 



y'" - (py)" + (qy)' -ry = o. 



The fact that %, z 2 , z 3 are solutions of (4) may be verified by substitution. 
The relation between (1) and (4) is entirely reciprocal in its nature, for 
equation (1) is also the adjoint to (4). Furthermore, the functions z lt z. 2) z 3 
are linearly independent solutions of (4) since otherwise the curve O would 
reduce to a point or straight line, and y lt y 2 , y 3 would be also linearly de- 
pendent. The reader is referred to Darboux, Thgorie des Surfaces, vol. II, 
book IV, chap. 5 for a discussion of the properties of the adjoint equation. 
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The identities 

2/1% + y&% + 2/3% = o, 

(5) y[z x + ?/ 2 z 2 + yfo, = 0, 

2/i'zi + 2/2% + 2/3% = 1» 

are obvious on substitution of the expressions given in (3) for z lf z 2 , z s . 
These equations serve also to define z u z 2 , z 3 in terms of y { , y 2 , y 3 . 

Returning now to the curve C, we see that it is a continuous curve in the 
projective plane (i.e. if the projective plane be mapped on the ordinary 
Euclidean plane so that P is a point of O in the finite plane, then C is con- 
tinuous at -P i n the ordinary sense) since y lf y 2 , y 3 are continuous functions, 
not all simultaneously zero ( W & 0) . Likewise O has a continuously turn- 
ing tangent in the projective plane (i.e. if the projective plane be mapped on 
the ordinary Euclidean plane so that P is a point of O in the finite plane, 
then the direction of the tangent line varies continuously along the curve at 
P„), since z x , ? 2 , z 3 are also continuous functions, not all simultaneously zero. 
Thus O is a regular curve in the projective plane. The fact that y u y it y 3 
have derivatives in x of all orders may be taken account of by saying that O 
is completely regular. 

The significance of the condition Wj± is that there are no points of in- 
flection. In fact, by definition, the tangent at such a point has contact with 
the curve of higher order than the first, and hence if 

atfOi + a 2 w 2 + a 3 w 3 = 

is the tangent line, one has simultaneously 

«i2/i + a 2 y 2 + a 3 y 3 = 0, 

»l2/l + a &2 + "32/3 = 0> 

«i2/i' + *&t + a s2/s = 0. 

This is possible if and only if W vanishes. The curve O may there/ore be any 
completely regular curve without any point of inflection. The reader will re- 
member that the only condition on y lt y. 2 , y 3 that they form linearly indepen- 
dent solutions of some equation (1) is the condition W & 0. 

If we consider the projective plane as lying in the ordinary Euclidean 
plane, the curve C possesses a tangent which rotates continuously in one 
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direction until the curve passes through the 'line at infinity' which is the im- 
age of any arbitrary line in the projective plane. It then returns with a tan- 
gent rotating in the reverse direction (consider for example a conic, which 
has no point of inflection). 
A transformation 

y = \y, (X*0), 

followed by a division by \, leaves (1) linear and of course does not change 
the curve O. By choosing 

X = e ~ U 1 "", 

the coefficient of the second term of (1) will reduce to zero and that equation 
will take on the simple form 

(6) y<" + qy' + ry = 0, 
to which the adjoint equation is 

(7) y«> + qy' + (q' - r)y = 0. 

It may be observed in passing that this last equation might also have 
been obtained from (4) by a like transformation in which however 

X = e V pdz 

A simplification introduced by giving our equations this normal form is that 
W reduces by (2) to a non-vanishing constant. 

From this point forth we shall often employ our equations in this simplified 
form. The curve O may then be looked upon as given either in normalized 
point coordinates by (6) or in normalized line co5rdinates by (7). 

The above interpretation of the solutions of an ordinary linear homogen- 
eous differential equation by means of a curve is a well-known one. The sec- 
ond and third chapters of Professor Wilezynski's Projective Differential Ge- 
ometry contain a full development of it. 

2. General Separation Theorem. The above geometrical repre- 
sentation furnishes us at once with the negative conclusion — not to be seen 
very readily without use of this representation perhaps — that it is possible to 
choose (a) the differential equation (1), and (b) three linearly independent 
solutions?/!, ?/2) y%> of the same, so that the zeros of y x , y 2 , y s will succeed each 
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other in any arbitrarily prescribed order. It will suffice to show that there 
exists a completely regular curve O without points of inflection which cuts the 
sides of the triangle of reference in any arbitrarily prescribed succession, since 
such a curve is the integral curve of some equation (1). The existence of such 
a curve may be readily seen as follows : Let G be drawn up to a point P not on 
a side of the triangle of reference ; the curve C may be extended beyond P 
without the introduction of points of inflection, so as to cut any assigned side 
of the triangle. The method of extension is given in the adjoining figure 




Fig. 1. 

(Fig. 1) . By a succession of such extensions, C may be made to cut w 1 = 0, 
w, 2 = 0, w s = in the desired order. 

It is only when we consider the solutions y x , y 2 , y 3 and z,, 2 2 , z s together 
that we obtain a Separation Theorem concerning the order of the zeros of 
these six functions. 

General Separation Theorem. Between any two successive zeros of 
t/i (or of Zi) there are an odd number of zeros of y k and z t together, where 
(i, k, I) is any permutation of (I, 2, 3). 

Proof. If we can prove that between any pair of successive zeros of y x 
there are an odd number of zeros of y. 2 and z 3 together, the theorem will follow 
by considerations of symmetry and duality. We may suppose without loss of 
generality that the line w l = of the projective plane goes over into the line 
at infinity in the Euclidean plane. The arc of C corresponding to values of 
x between the two values for which y, = is then a single branch of the curve. 
The line w. z = is some line in the finite plane. The geometrical interpreta- 
tion of the condition y 2 = is then that the curve O cuts this fixed line w 2 = ; 
the interpretation of the conditions = is that the tangent line is parallel to 
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the line w. z — 0, in other words, passes through the intersection of w 1 = and 



w>2 = 0. 



w a =0 




FIG. 2. 

Now (Fig. 2) if the arc enters on one side of w. t = in this plane and 
leaves on the other there are an odd number of zeros of y 2 and an even num- 
ber of zeros of z 3 ; otherwise there are an even number of zeros of y % , and an 
odd number of zeros of z s . In either case the theorem holds. 

Self-Adjoint Equations. If the integral curve O is self-dual, equations 
(6) and (7) are identical since the normal form (6) is unique. The con- 
dition for this is plainly , 

In this case, which possesses especial interest, z x , z 2 , z :i must be linearly 
expressible in terms of y u y 2 , y 3 : 

z \ = Cn2/i + c 12 y 2 + c 13 y 3 , 

z % - CnVi + Cjaya + c 23 ?/ 3 , 

z 3 = C 3\Vl + C 32?/2 + C 33y3- 

Since the solutions z lf z 2 , z 3 are linearly independent, the determinant of 
the coefficients is not zero. Hence one finds from the first equation (5), by 
substitution of the expressions above for z u z 2 , z 3 , that 

Cny! + c 21 y\ + c 3i yl + (c 2J + c 32 )y 2 y 3 + (c 31 + c Vi )y i y i + (c 12 + c^y^ = 0, 
the equation of a conic unless 

c ll — c 22 = c 33 = 0> c 23 — — c 32» c 31 = _ c 13< c 12 = — c 21- 

These conditions cannot all hold since, if they do, the determinant of the 
coefficients will be zero. The integral curve C must be an arc of the conic, or 



1911] DIFFERENTIAL EQUATIONS OF THE THIRD ORDER 111 

the entire conic. In the latter case as x increases from a to b, the correspond- 
ing point P may traverse the conic in one sense any number of times, making 
q complete circuits and a partial circuit. 

The above considerations load us at once to the following theorem : 
Separation Theorem. Self-Adjoint Case. There exist solutions of 
any equation (1) reducible to self-adjoint form which nowhere vanish in (a,b) . 
If i/i and y 2 are solutions of the equation with at least one zero in (a, b), the 
zeros of yi and y. 2 either separate each other singly or in pairs. 

Proof. In fact, there exist straight lines which do not cut the conic, and 
therefore solutions which do not vanish, since the condition 

a a/i + a-Mi + a a/3 * ° 

is equivalent to the condition that the conic does not cut the line 

a^Vx + a. 2 w. 2 + a 3 W 3 = 0. 

Also, of two straight lines which do cut the conic, say in M ly M 2 and iV" x , 
JV 2 , respectively, either Mi and M 2 separate JViand JV 2 on the conic, in which 
case the zeros of the corresponding solutions occur alternately as we traverse 
the conic ; or these pairs do not separate each other, in which case we have 
first two zeros of the one solution, then two of the other, and so on. 

If C is merely an arc of the conic no solution vanishes more than twice 
in (a, b) although the theorem still holds. 

Transformation of the Independent Variable. In obtaining the normal 
form (6) only a multiplicative transformation of the dependent variable, 

y = \y, (X^O), 

was made use of, while it is known that (1) also remains linear when any 
transformation of the independent variable 

x = 4>{x), (f#0), 

is made. From the standpoint of the geometrical interpretation, the second 
transformation changes the parameter of the curve. 

"What is the simplest form of the equation obtainable by a suitable com- 
bination of both transformations ? 

First, if the equation be reducible to self-adjoint form the locus of the 
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point (i/i, y 2 » !/i) i s a conic, and a suitable transformation of the variables 
will reduce y lt y. t , y 3 to the form 

y 1 = cos as, y 2 = sin x, y 3 - 1, 

no matter how many times the integral conic be traversed. Hence a normal 
form in this case is 

y>" + y' = 0, (a^xzk b). 

Likewise, since the equation (6) may be written 

y»' + qy'+ (L + fi^y = , 
one sees that the transformation which reduces the self-adjoint equation 

(8) ?/'" + W' + j?/ = 

to the preceding normal form will at the same time reduce (6) to an equation 

(9) y'" + y' + B y = 0, (a ^ x ^ b). 
The Forsyth-Laguerre normal form is 

y"' + liy = 0; 
in this case the self-adjoint equation (8) has been transformed to the equation 

?" = 
with solutions 1, x, x 1 . Since the equations 

W Y = 1, W 2 = X, ic s = X 2 

give us each point of the conic once and only once, this transformation is not 
proper unless the interval (a, b) is so small that the integral curve of (8) is 
an arc of the conic. 

3. Separation Theorem for Regular Intervals. The somewhat 
indefinite character of the first Separation Theorem is due in large measure to 
the fact that the integral curve C may have a very complicated form ; in par- 
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ticular it may have double points and doable tangents. A regular interval 
(a, b) is deiiiied to be an interval such that the corresponding integral curve O 
is either an oval, or has no double points or double tangents except where C 
touches but does not intersect itself. On account of the fact that double 
points and double tangents are dual notions the interval (a, b) will be regular 
for the adjoint equation (4) if it is regular for (1) itself, and conversely. It 
is obvious that any interval can be separated into a finite number of regular 
intervals. 

Separation Theorem for Regular Intervals. If (a, b) is a regular 
interval there will exist a solution ij\of (1) which is of one sign in (a, b) and 
a family of solutions c^y^ + c$y 3 such thai the zeros of any two members of the 
family separate each other. 

Proof. The theorem is clearly true if the integral curve is an oval, so 
that this possibility is at once disposed of. 

Assume first that the arc AB corresponding to the interval (a, b) does 
not cut some line of the projective plane, which is taken to be the line at in- 
finity in the Euclidean plane. On account of the hypothesis that (a, o) is 
regular, it is clear that AB will be spiraliform (Fig. 3). We are not 
excluding the possibility that A or B or both lie on the line at infinity. 




Fia. 3. 



We shall now show that there exists such a line in every case. Construct 
the tangent I at a point P of AB. As P moves from A to B along the arc 
AB, the line I will not at first cut AP. There must, however, be a least value 
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of x, say £, such that for slightly greater values the corresponding tangent 
I will cut AP; otherwise the tangent line at B does not cut AB and we are 
brought back to the previous case. Furthermore, the tangent corresponding 
to x = £ obviously passes through A. If this particular tangent line be chosen 
as the line at infinity, it is clear that AP will correspond to a single branch 
of AB, without inflection point or double point, of course, and that accordingly 
there exist lines I, which do not touch AP. Now take such a line I, as the line 
at infinity, and the arc AP will have the form indicated in Figure 4. 




Fig. i. 



From this figure it is clear that, since there are no double points of inter- 
section, the entire are AB must lie between the arc AP and the tangent at P. 
Hence this case is reduced to that first considered. 

Thus the integral curve O is spiraliform in all cases. It should be 
noted that an oval is to be looked upon as a limiting form of a spiral. 

Any line which does not have a point in common with AB will correspond 
to a solution y h which does not vanish for a & x ^ b. Likewise all the lines 
through a point in the interior of the integral curve (see shaded region in 
Fig. 4) will correspond to a linear family 

c 2 y 2 + c 3 2/ 3 

whose zeros separate each other. 

4. An Auxiliary Formula. We shall now derive an auxiliary 
formula upon which for the most part the further theorems will be based. 

Let y x , y 2 , y 3 be three linearly independent solutions of a given linear 
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differential equation (l),and z x , z. h z 3 the corresponding solutions of the adjoint 
equation (4), and let us write 

(10) </>(*, f) = yi(*)*i(f) + yi(*)«i(f) +y«(*)«i(f). 

The relations (5) are then 

(11) *(£,£) =*i(&£)=0. *i(f,f) = l » 
where our notation is 

^ l(a; * ^ = al ' M x ' *>■ = a?" • etc- 

Take now (1) in the special form (6) and consider the modified equation 

(12) y'" + qy' + (r + R)y = 0, 

which may be written in the form of a non-homogeneous linear differential 
equation 

( 13) y'" + qy' + ry =f= - By. 
The general solution of this equation may be written 

(14) y = c lVl + c$j t + c 3 y s + H>(x, £) /(f) df . 

In fact by direct differentiation and reduction by means of (11) one finds 

y' = ciy[ + <w* + c 3 yi + /&(*, S)A& d S> 

y" = c x y'{ + ctfS + c 3 yl> + jy. 2 {x, f )/f W* 

y'"= c iyi "+ c#H'+cirt'+/(z) + fax, f)/(f)rif. 

Substitute these expressions for y and its derivatives in (13). This equa- 
tion will then reduce to an identity by virtue of the fact that y x , y%, y s and 
(j>(x, f) are solutions of (6), which proves that the expression for y given by 
(14) is the general solution of (13). 
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The function 

V = Wi + 'hVt + c 0s 

is a solution of (6), and from the above equations we see that it satisfies the 
conditions 

(15) y(a) = V (a), y'(a) = ,'(«). y"(a) = V "(a). 

Thus we may rewrite (14) in the final form 

(i6) y = v -f}( x > £) R (t) y(f) *t> 

in which y and v are corresponding solutions of (12) and (6) respectively, 
related by equations (15). This is the desired formula, and it represents, 
of course, nothing more than the solution of a special non-homogeneous linear 
differential equation. 

5. Test for Regular Intervals. A regular interval (a, b) is said to 
be of the first kind if the tangent at A does not meet the arc .423 at anysecond 
point (see Fig. 4), and of the second kind if the tangent at B does not meet 
the arc AB at any second point. According to this definition a regular in- 
terval may be of both kinds. We shall take the interval (a, b) to be of the 
first kind, stating, however, the results for both cases. 

It follows at once from the definition, as has been noted, that if (a, b) is 
a regular interval for equation (1 ), it is also a regular interval for equation (4), 
the three solutions Zj, z % , z % of (4), given by (3), now being interpreted as the 
coSrdinates of a point. Let the corresponding point locus be the curve D. 
Then D is obtained from O by a polar reciprocation and a suitable projection. 
Furthermore, if (a, b) is regular and of the first kind for (1), no line through 
B is tangent to the curve C at any other point than B (see Fig. 4) . There- 
fore, in the curve D, no point on the tangent through B lies on the curve D, 
i. e. D is of the second kind. 

A necessary condition that (a, b) is a regular interval is that all the points 
of intersection of the tangent I at any point P of AB lie on one and the 
same side of P. Thus in the preceding figure (Fig. 4) all the points of inter- 
section of I lie in AP and none in PB. To prove the condition sufficient we 
have to show that, if it is satisfied, the integral curve AB can have no double 
points or double tangents except where the curve C touches but docs not inter- 
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sect itself. Suppose there were such a double point M = N where JVlies in 
MB. Then the tangent at M cuts the curve at JV", and the tangent at JV cuts 
the curve at M. This is impossible if the condition is satisfied. If there 
were such a double tangent at M and JV, where JV lies in MB, the tangent 
line in the vicinity of M would cut the integral curve in the vicinity of JV", 
and likewise the tangent line in the vicinity of JV would cut the integral 
curve in the vicinity of M. This is also impossible if the condition is sat- 
isfied. Thus the validity of the condition has been demonstrated. Further- 
more it is clear from the figure that (a, b) will be regular and of the first kind 
if no point of intersection of the tangent I lies in PB. 

The analytic phrasing of the geometric condition can easily be given in 
terms of the function <f>(x, £) defined in (10). For from the equations (11) 
it follows that the solution of (1) which corresponds to the tangent having a 
point of contact whose x is f, is §(x, f). Accordingly <f>(x, £) will not 
change sign for x > g, since the tangent at f does not cut the integral curve for 
points whoso x is greater than £. Moreover, <£(#, f), regarded as a func- 
tion of a;, has aminimum at the point x = £, by (11). Therefore the analytic 
form of the condition is the following : the interval (a, b) is regular and of 
the first (second) kind provided that 4>( x , £) is positive or zero for x> £(#<£). 

We are now in a position to develop the following test: If (a, b) is a 
reyular interval of the first (second) kind for the equation (ft) it is also a regular 
interval of the same kind for the modified equation (12), provided that the 
inequality It S (It S 0) obtains. 

Proof. The proof can immediately be made by means of the auxiliary 
formula (16) in which we take 

V = <f>(x, a), 

where a remains to be specified. This choice of n is possible since <f>(x, a) 
is a solution of (6) . We obtain then 

(17) y = $(x, a) - j%(*, ?) i?(£) y (*)<**, 

when y is that solution of the modified equation (12) which by (15) satisfies 
the conditions 

(18) y(a) = <f>(a, a) = 0, y'(a) = ^(a, a) = 0, y"(a) = ft(a, a) = 1. 
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Accordingly y is precisely the solution of the modified equation (12) which 
plays the same role for (12) that <f>(x, a) does for the equation (6). There- 
fore if (a, b) is not regular and of the first kind for the modified equation 
there exists some a such that the corresponding y changes sign for x > a. Let 
x be the least value of x for which this is true. Then we have the equation 

o = <f>(x , a) - jT?(xo, n R{Z) y(f ) *£• 

But by hypothesis (a, b) is regular and of the first kind so that <f>(x , a) 
and (f>(x , l;)(x > f) are positive or zero. Moreover i?(f) is negative by 
hypothesis, and y(f) is positive for f between a and x by (18). Thus we 
have the sum of two positive or zero quantities equal to zero. Of these two 
quantities the second given by the definite integral must be negative ; for the 
first factor <f>(x , f) of the integrand can only be zero for isolated values 
of x since <j>(Xo, !• ) is a not identically zero solution in f of the adjoint equation 
to (6), given by 

yiO«o)zi(£) + yifoHCl) + ys(^o)^(l) ; 

likewise the third factor of the integral cannot vanish except at x and 
£ on account of the way in which x was chosen ; finally the remaining factor 
R(%) may be taken to be not identically zero. Hence the above equation is 
impossible. 

The following special case of the test is important: The interval (a, b) 
is a regular interval of the first (second} kind for an equation 

(19) y'" + gy' + (I + R) y = 

provided that R ^ 0(22 & 0) in the interval (a, b). 

In this case it is only necessary to identify the equation ((5) with the self- 
adjoint equation : 

</" + gy'+f'y = o, 

for which any interval (a, b) is regular and of both kinds, and to state the 
test for this special case. 
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Since any equation (6) can be written in the form (19) we have a method 
of separating any interval into a sequence of regular intervals the first of one 
kind, the next of the other, and so on, the end points of these intervals being 
given by those roots of R = at which i2 changes sign. 

The regular intervals thus chosen have a special geometric property which 
however I shall not stop to prove : The successive osculating conies to the in- 
tegral curve lie one within the other at all points of such an interval, and at the 
end points the curve O and the osculating conic hyper osculate. 

6. Intervals of Oscillation. Comparison Theorem. It remains 
to develop the nature of the dependence of solution and coefficients. Anal- 
ogous to the notion of the interval of oscillation as the distance between suc- 
cessive zeros of a solution of an equation of the second order, is the notion of 
forward and backward intervals of oscillation. 

Consider the totality of those solutions of (1) which vanish at x = a j 
these solutions will form a two parameter linear family. Suppose that they all 
vanish again in the interval (a, b). The least interval (a, /3), where /3 is 
greater than a, which has the property that all these solutions vanishing aj 
x = a will vanish again in (a, /3) is called the forward interval of oscillation 
at x = a. Likewise the least interval (/3, a) where is less than a, which has 
the property that all these solutions vanish again in (/3, a) is called the back- 
ward interval of oscillation at x = a. 




Fig. 5. 

It may happen that not all of the solutions which, vanish for x — a vanish 
again for x > a ; in this case a forward interval of oscillation at x — a does not 
exist. Likewise a backward interval of oscillation may not exist. 

Let us now consider a regular interval (a, b) which, as before, is taken to 
be of the first kind. The integral curve AB will then be a sphaliform curve 
with the point B in its interior (see Fig. 5). For such an interval (a, b) a 
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forward interval of oscillation never exists, unless AB is an oval or touches 
itself. Thus in the figure let P be any point whose x is a ; then the line PL 
corresponds to a solution which does not vanish for x > a, since PL does not 
meet PB at any second point. 

The backward interval of oscillation may however exist. Thus in the 
figure, the tangent to the curve AB at P first meets the arc AP&t Q ; if the 
point Q corresponds to x = ft, then (ft, a) is the backward interval of os- 
cillation at a. Any straight line through P meets QP again and this is not 
true of any shorter arc PQ'. Since <$>(x,a) corresponds to the tangent at a, 
the quantity ft is the greatest root of <fr(x,a) = less than a. If the tangent at 
Pdoes not meet AP again, the backward interval of oscillation does not exist. 

It is apparent from this figure that any solution vanishes at most in two 
distinct points in a backward interval of oscillation (ft, a) on such an interval 
(a, b). Furthermore, if a given solution y vanishes for any value of x > a, 
it will necessarily vanish in (ft, a). This is obvious from the figure also. 

Similar facts are true of forward intervals of oscillation if (a, b) is of the 
second kind. 

Comparison Theorem. If (a, b) be a regular interval of the first 
{second) kind for the equation (6), and )/5S0(iJ^0), the backward 
(forward) interval of oscillation (&,<*) [(a,/8 x )] at x = a for (12) is smaller 
than the like interval (ft, a) [(«,#)] for (6). 

Proof. To prove this theorem we shall again employ the auxiliary form- 
ula (16). As in a preceding proof we shall take 

t) = 4> (x, a) 

and thus obtain again (17). Suppose now if possible that the theorem is not 
true, so that we have /3 S ft. As we have already noted, the function analo- 
gous to (f)(x, a) for the modified equation (12) is precisely y of (17), which 
must then vanish for x = & but not for X < x < a, by the analytic condition 
obtained above, since (a, b) is regular for (12) also. The two functions y 
and (f>(x, a) will therefore be positive throughout (ft, a) since neither vanishes 
within (ft, a) and since both are positive in the vicinity of a by (1 1). Put 
now x = ft in equation (17). It reduces to 

y(ft) =J%(ft, £) JJ(£) y(£)rff, 
where <j>(ft, f) is positive. But i?(f) is negative or zero. Therefore in this 
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relation are equated quantities, one positive or zero, the other actually nega- 
tive. This is impossible and the theorem must hold. 

We note the following result contained in the above theorem : The back- 
ward {forward) interval of oscillation for the equation 



y'" + qy' + (^+ R)y = o, 



where R ^ (R S 0) is smaller than that for the self -adjoint equation obtained 
by taking R = 0. 

ZU- 




FIG. 6. 

Any equation (1) may be reduced to the normal form above by a known 
transformation. In the present section we have developed the dependence of the 
solutions on R, by showing that if R is negative a decrease in R decreases the 
backward interval of oscillation. In the figure (Fig. 6) we have repre- 
sented by a full line the integral curve O of this equation and by a broken line 
the integral curve D of the modified equation. If (/3, a) represents the backward 
interval of oscillation at a for the curve C and (&, a) the similar interval 
for D, we have /8 X > $. For the sake of the comparison the curves 
have been taken with the same initial point and direction. It would 
now be desirable to give a comparison theorem which stated, in simple terms, 
the dependence of the interval of oscillation on the function q. If we increase 
q, however, it may be shown that the interval of oscillation may be either in- 
creased or diminished in length. In the self-adjoint case the dependence is 
more simple : 

Comparison Theorem. Self-adjoint case. Given two equations reduci- 
ble to self-adjoint form, in particular to 

y'" + w' + -f y = o, 

(20) 

y'" + q*y' + q fy = o, 
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such that qi > q 2 in (a, b) ; then the backward (forward) interval 0/ oscillation 
of the first of these equations is less than that of the second. 

Proof In order to prove this theorem we shall prove first that the back- 
ward interval (/S, a) of oscillation of any self-adjoint equation 



(21) 2/ 



'// 



+ qy + -^r 1 v = ° 



is precisely the interval of oscillation (fi,a) for the linear differential equation 
of the second order 

(22) z" + I z = 0. 

Let y be the solution of the equation (21) such that 

(23) y(a)=y'(a)=0, y"(a) = 1. 

That is, y is chosen so that it is equal to <j>(x, a) by (11). As has been shown, 
the backward interval (#, a) of oscillation at x = a is given by the largest 
root of y = which is less than a. 

Now we have clearly, by triple integration by parts, 

/ yy'"dx = yy" - y' 2 + y"y - i y'"ydx, 

and by a single integration by parts, 

y(qy') dx = qy* - / {qy' + q'y)ydx. 

Therefore we have 

fy(y"' + qy' +\ y)dx = 2yy" - y'* + qy 2 - J \y'" + qy' +|' y)y dx. 

Kemembering now that y is a solution of equation (21) which furthermore 
satisfies conditions (23) we obtain 

2yy" - y<* + qy 2 = 0. 
Make the substitution 

y = z\ y' = 2 zz\ y" = 2 (zz" + z' 2 ) , 
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in this last equation, and, on the removal of a factor 4s 3 it will reduce to the 
linear differential equation (22). It is clear that, as thus determined, z will 
be a solution which vanishes at x = a. And # will be determined of course as 
the greatest value of x less than a for which z = 0. Hence the backward 
{forward) interval of oscillation at a of the given self-adjoint equation (21) 
is given as an interval between two successive zeros, of a solution of the asso- 
ciated linear differential equation (22) of the second order. 

The theorem we wish to prove is now an obvious consequence of the 
Comparison Theorem stated in the introduction. 

The above reduction is not new. 

7. Oscillatory and Non-Oscillatory Solutions. The preceding 
paragraphs afford practical methods for obtaining an idea as to the distribution 
of the zeros of a solution of (1) provided the solution is known to vanish. 
The following theorem often enables one to determine whether a given solution 
is oscillatory or non-oscillatory in the interval (a, b) i. e. whether the solution 
vanishes or not. 

Theorem. Let (a, b) be a regular interval of the first (second) kind for 
the equation (6), and suppose that It SO (R § 0). Let y and y be solu- 
tions of (6) and (12) respectively, so related that 

y{a) = v (a), y'(a) = y' (a), y"(a) = y"(a). 
The inequality 

(24) \y\ < \y\ 

obtains for values of x less (greater) than a and within the backward 
(forward) interval of oscillation at x = a until a value ofx is reached for which 
y = ; the inequality | y | > | r\ | obtains for all values of x greater (less) than a 
until a value of x is reached for which y = 0. 

Proof If the inequality (24) fails to hold force < a, it will certainly 
fail to hold before an x is reached for which one has y = 0. In fact, otherwise, 
for this x we have y = 0, by (24), and the theorem does not say anything 
about y for lesser values of x. Therefore if the first part of the theorem does 
not hold, the inequality 

(25) \y\ S \y\ 

is true for some x < a, say cc , in the backward interval of oscillation (/3,a), 
such that neither y nor y vanish in (x , a) . 
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But y and y have precisely the properties that they had in the auxiliary 
formula (16) so that 

y(x ) _ V (x ) = [<]>(x , £) Rtf) y(l) df. 

The left h ind member in this equality is of the same sign as y(x ) or zero, since 
(25) holds at x = x . The right hand member is of opposite sign to y(x ) , since 
R(l-) is negative and <f>(x , £) is positive, x and f being contained within the 
interval of oscillation (#, a) (see Fig. 5). This equality cannot hold and the 
first part of the theorem is true. 

The second part of the theorem is proved in an analogous way. An 
equation 

y(x ) - V (x ) =-J%(x , £) JB(f) y(£) d £ 

will be obtained if the theorem is not true, in which x > a is an x for which 
|j/|^|ij| while neither y nor?; change sign in (a, a: ). In this case too R(%) 
is negative and<|>(a; , f) is positive. As before, the two members of this equal- 
ity would have opposite signs. 

The relation of the two solutions is given in the figure (Fig. 7). 

1Z 




Fro. 7. 



The following application is at once obvious : If y and y are solutions 

of the equations y'" + qy' + ^- y = and of ?/"' + qy' + (~ + R) y = 

respectively, where R ^ (R § 0), and are so related that 

y(a) = y (a), y'(a) = V '(a), y"(a) = y"(a), 
then \y\ < \y\ forx < a(x > a) and \y\>\y\ for x > a(x < a) until a value 
of x is reached for which y = and y = respectively. 
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In fact the interval (a, b) under consideration is regular and of either 
kind for the self-adjoint equation. 

8. An Example. It is possible in many cases to form a good idea of 
the character of the solutions by means of the preceding results. As an illus- 
tration we consider the equation 

y'" + if + xy = 0, 
to which the adjoint equation is 

z"> + z < — xz = 0. 

We choose y lt y 2 , y s as the principal solutions at x = : 

yi(0) = 0, y[ (0) = 0, y{> (0) = 1, 

y 8 (0) - 0, y> % (0) = 1, y' 2 (0) = 0, 

7/3 (0) =1, y{ (0) = 0, y'i (0) = 0. 

The adjoint solutions (see equations (3) then satisfy the conditions 

*(0) = i, «I(0) = o, 4'(0)=-i, 
z% (0) = o, 4 (0) = -i, 4 (0) = o, 

z 3 (0) = 0, 4 (0) = 0, 4 (0) = 1. 

The interval x < is regular and of the first kind (test for regular in- 
tervals, special case) and the interval x > is regular and of the second kind. 
Each of the corresponding parts of the integral curve is therefore spiraliform 
with the point A corresponding to x = in its interior. 

Furthermore it is clear that y u y 3 , z u z 3 , are even functions of x and that 
y 2 , n are odd functions, since the equation is unchanged if —x be substituted 
for x. The integral curve must be symmetrical with respect to w 2 = 0. Hence 
the integral curve has the form indicated in the schematic figure (Fig. 8). 

The zeros of y l and y 2 and indeed of any two members of the linear 
family c x y x -f- c 2 y.^ will alternate for x > and for x < (Separation Theo- 
rem for Regular Intervals) . 

In order to consider the question of the distribution of the zeros of y u 
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■i/i for x > we note that the forward interval of oscillation will decrease with 
increase of x and always be less than that of 

y'" + y = 

which is precisely 2ir (Comparison Theorem). Hence as x increases the 
successive zeros of y x and of y % come at intervals none exceeding 2ir and 
tending to zero. 

The function y 3 also vanishes for x = « x > 0, as a direct computation 
shows. The successive zeros of y A {x) will fall of course at intervals less 
than 2-ir which tend to zero as x increases. 



,Wf=0 




FIG. 8. 

9. The Theorems Of Liouville. The only known theorems of the 
kind we have been considering concerning the solution of the third (and 
higher) order equations appear to be those due to Liouville.* 

The type of the equation considered is 

d __ d T d 

ix- K rx L dx^ + ^ = ' 

where K, L are positive in (a, b) and \ is a parameter. Now when X = 
one finds by direct integration that 



<KX, a) = j'j. p^dx. 



*Liouville's Journal, vol. 3, pp. 561-614 (1838). 
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This function does not vanish for x ^ a. Hence for X = the interval (a, b) 
is of the first and second kind, the corresponding arc AB has the shape of 
a part of a convex oval. Suppose now that the equation is reduced to the 
form (6) and that X is made to increase. By the test for regular intervals and 
the Comparison Theorem, the interval (a, b) remains regular of the second 
kind and the forward intervals of oscillation decrease. In addition, we know 
by the theorem of §7 that if we fix y, ?/, y" at x = a, the first greater value 
of x for which y = will decrease with increase of \. 

The results of Liouville are still more definite, such as the following : If 
we consider a particular solution y(x) for which at x = a 

y=a > Tx y ~ h > dx L d&y = c ' a>0 > 6>0 ' c>0 ' 

all of the greater values of x for which y vanishes must decrease with increase 
of X. 

This result is equivalent to his oscillation theorem which states that there 
exists one and only one value of X, say X , for which y vanishes at b and not 
within (a, b) , and one and but one value X lf such that y vanishes at b and once 
within (a, b) and so on. 

Liouville's results belong to another category than those which we have 
obtained. We have given only theorems concerned with any arbitrary in- 
crease of the coefficients. A special kind of variation was considered by 
Liouville which is however of great importance. 

Princeton University, 
June, 1910. 



